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Qh' We calculate the heavy quarkonium production in double pomeron exchange 

' processes in perturbative QCD by using two-gluon exchange model. For the 

On . P-wave xj productions, we find xi ^i^d X2 production amplitudes which van- 

^ \ ish in the forward scattering limit. We also calculate direct J/ipiT) + 7 

I production in the same approach, and these direct contributions are much 
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smaller than the feeddown contributions from the P-wave states. 
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In recent years, there has been a renaissance of interest in diffractive scattering. These 
diffractive processes are described by the Regge theory in terms of the Pomeron (IP) exchange 
|1|]. The Pomeron carries quantum numbers of the vacuum, so it is a colorless entity in QCD 
language, which may lead to the "rapidity gap" events in experiments. However, the nature 
of Pomeron and its reaction with hadrons remain a mystery. For a long time it had been 
understood that the dynamics of the "soft pomeron" is deeply tied to confinement 0. 
However, it has been realized now that much can be learned about QCD from the wide 
variety of small-x and hard diffractive processes, which are now under study experimentally. 
In Refs. [0,^, the diffractive J/il) and T production cross section have been formulated in 
photoproduction processes and in DIS processes in perturbative QCD. In the framework 
of perturbative QCD the Pomeron is represented by a pair of gluon in the color-singlet 
state. This two-gluon exchange model can successfully describe the experimental results 
from HERA §. 

This two-gluon exchange model has been extended to the hadron hadron collisions for 
the single diffractive processes 0J^, where the so-called coherent diffraction mechanism 
plays an important role. In these processes, the Pomeron represented by a color-singlet 
two-gluon system is emitted from one hadron and interacts with another hadron in a hard 
process, in which the two gluons are both involved. These coherent diffraction processes 
are very interesting for the study of hard diffraction mechanism in hadron collisions, where 
there exist nonfactorization effects in the leading twist 0. 

In this paper, we will further extend the above ideas to the double Pomeron exchange 
(DPE) processes. In particular, we will calculate heavy quarkonium productions in DPE 
processes at hadron colliders. Apart from the single diffractive processes, the DPE processes 
at pp or pp colliders are also useful to the studies of the diffraction mechanism and exploration 
of the nature of the Pomeron. In the DPE, both incoming hadrons are quasielastically 
scattered and therefore emerge at very large Feynman xf, |a;F,i|, \xf,2\ > 0.9. And the final 
states system in the center region is produced from the Pomeron-Pomeron interaction. Due 
to the color-singlet (Pomeron) exchange, the final states system can be separated from the 
two beam particles (proton or antiproton) by large rapidity gaps. These DPE processes are 
quite attractive because it has been long time proposed to search Higgs meson in double 
rapidity gap events at proton colliders [|10|-[T6[|. An obvious advantage of the rapidity gap 
approach is the spectacularly clean experimental signature . 

In the two-gluon exchange model, Xc and xt states productions in DPE processes have 



been estimated in Ref. |T^, where they approximately take the gluons coupling to the x 
states in gg xj from their decay widths in xj ~^ 99 processes. In this paper, we will 
directly calculate these amplitudes in terms of the wave functions derivatives at the origin 
of the P-wave heavy quarkonium bound states in the framework of color-singlet model by 
using the nonrelativistic approximation, and we also take into account the strong correla- 
tions between the two incident gluons polarizations, which will enter into the production 
amplitudes explictly. As a result, we find vanishing amplitudes for xi and X2 states in the 
forward scattering limit. We will also calculate J/il) + '^ production in DPE processes in the 
same approach. 

The double pomeron exchange processes x production in 

p + p{p) -^p + Xj + PiP) (1) 
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is plotted in Fig. 1 in the two gluon exchange model, where the hard subprocess gg ^ M is 
initiated by gluon-gluon fusion and the second t-channel gluon (with transverse momentum 
Qt) is needed to screen the color flow across the rapidity gap intervals. In our calculations, 
we express the formulas in terms of the Sudakov variables. That is, every four-momenta ki 
are decomposed as, 

ki = aipi + (3iP2 + kiT, (2) 

where pi and p2 are the momenta of the incident hadrons (proton or antiproton). For high 
energy process, we will have pi = 0, P2 = 0; 2pi ■ P2 = s. ai and are the momentum 
fractions of pi and p2 respectively, kix is the transverse momentum, which satisfles 

kiT ■ Pi = 0, kiT ■P2 = 0. (3) 

For the momenta of the incident gluons qi and q2, we have the following decompositions, 

qi = xipi + giT, q2 = X2P2 + q2T- (4) 
And then, the momentum of the produced x state will have the following form, 

P{x) = qi + q2= XiPi + X2P2 + Pt, (5) 

where the x state transverse momentum Pt = qir + q2T- In the following, we will calculate 
the amplitudes in the forward scattering limit. That is to say, the momenta transfer squared 
between the flnal state hadrons and initiate hadrons are set to be zero, ti = {pi — p'lY = 
and t2 = {p2 — p'2)^ = 0. Under this limit, we will have glr = —q2T = Qt, and then 
Pt{x) = 0- With the above kinematical considerations, we can formulate the forward 
scattering amplitudes for — > p + x + p as follows, 

Mixj) = — j :^^^fg{x^,x[,QT,M)fg{x2,x'„QT,M) x Aj, (6) 

where fg are the unintegrated off-forward (skewed) gluon distribution functions. They can 
be related to the conventional gluon density as [|18[ 

/,(x, x', Qt, M) = R^^^l^jTiQr, M)xg{x, Ql)]. (7) 

The factor Rg takes into account the skewed effects of the off-forward gluon density compared 
with the conventional gluon density in the region of x' -C x. The bremsstrahlung survival 
probability can be found in [l^, which will reduce to the conventional Sudakov form 



factors in the double logarithmic limit. 

The expressions Aj in the amplitude (H) include the xj bound state effects (the wave 
functions), and will also depend on the polarizations of xj states. For J = 0, we have the 
following simple result. 



l2MaM 1 ... 
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where R'p{0) is the radial wave function derivative at the origin for the P-wave x state of 
charmonium or bottomonium. For J = 1, 



8 X 1287ra, 



9v^s(M2 + QlY V 47rM 



R'p{0)e 



Qtp1p2 



(9) 



where e(j=i) is the polarization vector of Xi state. From the above expression we can see 
that Ai oc Q^, so the amplitude for J = 1 state will be zero after angular integration of @ 
over Qt- That is to say, there is no contribution for Xi production in the DPE processes 

(S). 

For J = 2, 



Here 6(^=2) polarization tensor for X2 state, which has the following properties. 



For the angular integration of (BI) over Qt, we have the following identity. 



d'QTQf^Q^T = ^ I dQlQlg%\ 



(10) 



(11) 
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where the transverse part of g^y is given by 



So, after integrating the azimuth angle of Qt, the expression in the brackets 
will become 



of Eq. (Iig 



'^dfiuQ'r) 



(13) 



With the above results, we find that the amplitude for X2 will be equal to zero because of 
the identities of Eq. (pA]). 

The vanishing of Xi ^^id X2 production amplitudes are direct consequences of the for- 
ward scattering limits in the DPE processes formulated in Eq. (|^). Under this limit, the 
polarization vectors of the two incident gluons are strongly correlated, i.e., Cj oc qiT with 

QlT 

On the other hand, we know 



-g*2T |19-21,18|. Because of these strong correlations, the final state {xj state here) 
in the hard DPE processes must be in the 



state 21 



that the state of = decouples with two real gluons in the nonrelativistic limit for the 
heavy quarkonium system |2^. By explicit calculations in the above, we show that the DPE 
amplitudes for xi and X2 productions have the properties in (|[) and (p!OD, and finally give 
vanishing contribution to these states. We note that these properties are common results in 
the two-gluon exchange Pomeron model (either perturbative or nonperturbative two-gluon 
exchange models). However, in the other Pomeron model, such as D-L model ||23|, the pro- 
duction amplitudes for xi and X2 do not vanish p4| , ^ . This is because in this model the 
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Pomeron couples approximately like a C = +1 photon, and the DPE xj amplitudes were 
adopted from the analogous process of xj ~^ 1*1*-, which will give the same magnitude 
production rates for xi and X2 as that of Xo p4| , p5| . 

Finally, the differential cross section of Xo production will be 



da 



\Ml_ 

16^77'^ 



e^^^e^^^dtidt2dy, 



(14) 



where y is the rapidity of xo state, ti is the momentum transfer squared at the proton 
(antiproton) vertices, and h is slope for the proton form factor, which will be taken as 
h = 4:.0GeV~^ in the following numerical calculations. After integrating ti and t2, we get 



da 

dy 



dQl fg{xu x[, Qt, M)fg{x2, x^, Qt, M) 



n 2 



my 



(M2 + 2Ql 
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(15) 



In Fig. 2, we plot the differential cross sections ([ISD for xca and Xfeo(l-P) at the Fermilab 
Tevatron RUN II, where we adopt the rapidity gap survival probability to be 5*^ = 0.05 as 
in Ref. |T^. For the masses and wave functions we use 



M(xco) 
M(xw) 



3.4GeV, \B!, 
-- 10.2GeV, 



km 



0.075GeV^, 
= lA2GeV\ 



For the gluon distribution function, we use the GRV LO parameterization After inte- 
grating over the rapidity, we get the following cross sections for double- diffract ive productions 

of chico and Xbo as 



o"(Xco) = 735 nb, a{xbo) = 0.88 nb 



(16) 



. Comparing our results with those of 
with 



1^ 



we find that our result for Xco is roughly agreeing 

This is because we use 



However, our result for Xbo is much larger than that in 
the potential model calculations for the wave functions, while they use lattice calculations 
for the decay widths for xt states. So, for bb bound states, there exist large theoretical 
uncertainties associated with the bound states properties. 

For the experimental observation, with the high resolution missing-mass measurements 
at the Tevatron (AM ^ 250Mey |T^), we can identify all of the x states in the DPE 
processes. Furthermore, we can also measure xo productions in DPE by detecting their 
radiative decays to J/tp (T) plus photon. The branching ratio of xo ~^ J I ''Pi is about 
6.6 X 10"'^, so the final cross section for J/ip production from Xco feeddown decays is about 
c"(Xco) ^ Br{xco J/i^l) X Br{J/ip /^"""/i") ~ 300 pb. Having this large cross section, 
J/ip production in the DPE processes can be well measured with large number of lepton pair 
(yU+z/^ or e~^e~) events. For P-wave bottomonium states, we know that the branching ratio 
of Xw(l-P) ^ '^{^S)-y) is very small (the upper bound is BR < 6%), so it may be difficult 
to detect T(15') from Xboi^P) decays. However, for 2P state Xfeo(2-P), we have not so small 
branching ratios decaying into T(15') and T(25'). So, from the combining branching ratios, 
we will have 

a{xbo{2P)) X [Br{xbo{2P) ^ T{2S)j ^ + Br{xbo{2P) T(1^)X ^ /iV" 

^ 0.8 pb. 
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In the above estimate, we assumed that the production cross section of Xboi'^P) is the same 
as that of Xboi^P) because their wave functions at the origin are roughly the same ||2^. So, 
T states from the feeddown decays of Xfeo(2-P) can also be observable in the DPE processes 
at the Tevatron II with integrated luminosity of 15fb~^. 

Apart from the feeddown contributions from xj states, J/ if) and T can also be produced 
directly from the DPE processes in 

p + p{p) ^ p + J + p{p)-, (17) 

for which we plot in Fig. 3 the typical diagram. For the final produced J/z/', we will have 
the Sudakov variables as 

P = a^pi + (3^p2 + Pt, 

and for the outgoing photon, 

p^ = aspi + P3P2 + Pt- 

In the forward scattering limit, because qiT + q2T = 0, the transverse momentum of the 
photon will be equal to the transverse momentum of J/tp in balance, i.e., pj' + Pt = 0. The 
onshell conditions for the outgoing J/ip and photon will give the following relations for the 
associated Sudakov variables, 

a^P^ = — = — , asPs = — . (18) 

s s s 

From the momenta conservation, we also have 

Using the above Sudakov variables, we derive the differential cross section for J/ip'j 
production in the DPE processes. 



da \M\ 



dy^dy^dPPT 16%652 • 



(19) 



where and y^ are the rapidities of final state J/ip and photon, and and as can be 
related to rapidities, 

"^T „, Pt y 

The amplitude for J /ip + production is more complicated than those for xj productions. 
In the limit of Qt — > 0, we can simplify the amplitude and get the amplitude squared as 

2 X lQ^7:^ala,elalalm\Rsm\ 2 ...^ 
' ' ?>xl{asml + a^P^Y{mlY ^' ^ ' 

where 

= j ^^/,(xi,a;;,QT,M)/,(x2,x^,QT,M). (21) 
So the cross section will be 
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da _ 2TT^a^^aeelalalm\RsiO)\\^^ ^^^^ 



2 \2 9' 

t) 



In Fig. 4, we plot the differential cross sections da/dydPT as functions of Pt a.t y = 
for J/ip and T(IS') at the Tevatron RUN II. For the numerical inputs for masses and wave 
functions, we use p6[] 

M{J/tlj) = 3.1GeV, \Rs{0)\l = OMGeV^, 
M(T(15)) = 9A6GeV, \Rs{0)\l = 6.48^6 V^^ 

From Fig. 4, we see that the cross section for J/ip drops very rapidly with Pt, and the cross 
section for T drops much more slowly than that of J/ip. At low Pt, the cross section of J/i/j 
is much larger than that of T, but at large enough P^ it will be smaller than T. However, 
comparing the results in Fig. 4 to those in Fig. 2, we find that the cross sections for S'-wave 
quarkonia are much smaller than those for P-wave states, even after the latter cross sections 
have been multiplied by the decay branching ratios of x J/i^C^)^- So, we conclude that 
the dominant contributions to S'-wave quarkonium production in the DPE processes are 
the feeddown contributions from P-wave decays. This is quite different from the inclusive 
production processes or the DPE J/ip production in the IS Pomeron model PB|, where the 
direct production is the dominant source to the prompt J/ip productions. 

In conclusion, in this paper we have investigated the heavy quarkonium productions in 
the DPE processes at hadron colliders in perturbative QCD. We found that in the forward 
scattering limit, the amplitudes for xi and X2 vanish due to the strong correlation between 
the polarization vectors of the two incident gluons. We also calculated the direct S'-wave 
quarkonium production in the DPE processes associated with a photon, and found that the 
direct contributions are much smaller than the feeddown contributions from the P-wave 
states. 
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Figure Captions 



Fig.l. Sketch diagram of the DPE processes xj production at hadron coUiders: (a) Pomeron- 
Pomeron fusion into xj; (b) Double Pomeron exchange in the perturbative two-gluon ex- 
change picture. 

Fig. 2. The differential cross sections da/y for Xca and Xho in DPE processes at the Fermilab 
Tevatron RUN II. 

Fig. 3. A typical diagram for the direct J/il> (T) production in the DPE processes. 

Fig.4. The differential cross sections da/dydpr for J/ip and T production in the DPE at 

the Tevatron RUN II. 
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(a) 




Fig.l Sketch diagram of the DPE processes xj production at hadron col- 
liders: (a) Pomeron-Pomeron fusion into xj\ (b) Double Pomeron exchange 
in the perturbative two-gluon exchange picture. 




Fig. 3 A typical diagram for the direct J /ip (T) production in the DPE pro- 
cesses. 



